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Abstract 

We use an analytical model to describe the magnetocrystalline anisotropy energy (MAE) in solids as a function of band 
filling. The MAE is evaluated in second-order perturbation theory, which makes it possible to decompose the MAE into a sum 
of transitions between occupied and unoccupied pairs. The model enables us to characterize the MAE as a sum of contributions 
from different, often competing terms. The nitridometalates Li 2 [(Lii_ x T x )N], with T=Mn, Fe, Co, Ni, provide a system where 
the model is very effective because atomic like orbital characters are preserved and the decomposition is fairly clean. Model 
results are also compared against MAE evaluated directly from first-principles calculations for this system. Good qualitative 
agreement is found. 


I. INTRODUCTION 

Magnetocrystalline anisotropy is a particularly impor¬ 
tant intrinsic magnetic property pQ. Materials with per¬ 
pendicular magnetic anisotropy are used in an enormous 
variety of applications, including permanent magnets, 
magnetic random access memory, magnetic storage de¬ 
vices, and other spintronics applications. [2H5] 

Modern band theory methods have been widely used 
to investigate the magnetocrystalline anisotropy energy 
(MAE) in many systems [6j [7]. The MAE in a uni¬ 
axial system can be obtained by calculating the total- 
energy difference between different spin orientations (out 
of plane and in plane). However, MAE is usually a small 
quantity and a reliable ab initio calculation requires very 
precise, extensive calculations. Moreover, MAE is, in 
general, harder to interpret from the electronic structure 
than other properties, such as the magnetization. MAE 
often depends on very delicate details of the electronic 
structure |S]. Using perturbation theory, the MAE can be 
decomposed into virtual transitions between different or¬ 
bital pairs. In practice, the d bandwidth is large enough 
that it is nontrivial to meaningfully resolve the MAE into 
orbital components and predict its dependence on band 
filling. 

The magnetocrystalline anisotropy originates from 
spin-orbit coupling (SOC) [2] or, more precisely, the 
change in SOC as the spin-quantization axis rotates. In¬ 
cluding the relativistic corrections to the Hamiltonian 
lowers the system energy and breaks the rotational in¬ 
variance with respect to the spin-quantization axis. Here 
we refer to the additional energy due to the relativistic 
correction as SOC energy or relativistic energy E r . MAE 
is a result of the interplay between SOC and the crystal 
field [101. The MAE and change in orbital moment on 
rotation of the spin-quantization axis are closely related. 
We describe this below and denote them as K and Kl, 
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respectively. Without the SOC, the orbital moment is 
totally quenched by the crystal field in solids. Except for 
very heavy elements such as the actinides, SOC usually 
alleviates only a small part of the quenching and induces 
a small orbital moment relative to the spin moment. For 
3d transition metals, SOC is often much smaller than 
the bandwidth and crystal field splitting, and thus can 
be neglected in a first approximation. While the E r is 
generally small, its anisotropy with respect to spin rota¬ 
tion is often even orders of magnitude smaller. 

Recently, it had been found that a very high magnetic 
anisotropy can be obtained in 3d systems such as lithium 
nitridoferrate Lio [(Lii _^FeAN] pTHT4] , which can be 
viewed as an a-L^N crystal with Fe impurities. As 
found both in experiments [15] and calculations HHIE] us- 
ing density functional theory (DFT), the L^Lii-^Fe^N 
system possesses an extraordinary uniaxial anisotropy 
that originates from Fe impurities. The linear geom¬ 
etry of Fe-impurity sites results in an atomic like or¬ 
bital and then a large MAE. As found in both x-ray 
absorption spectroscopy pT and DFT calculations 11 
na, 3d ions T have an unusually low oxidation state 
(+1 ) in Li 2 (Lii_ x T :c )N for T= Mn, Fe, Co, and 
Ni. Recently, Jesche et al.[ Id developed a single-crystal 
growth technique for these systems and directly ob¬ 
served that the MAE oscillates when progressing from 
T=Mn—)-Fe—)>Co—)*Ni. [ITS] Electronic structure calcula¬ 
tions also show that the atomic like orbital features 
are preserved for different T elements. Considering the 
rather large MAE and well-separated density of states 
(DOS) peaks in this system, it provides us with a unique 
platform to investigate the MAE as a function of band 
filling. 

Li and N are very light elements with s and p electrons, 
respectively. They barely contribute to the MAE in 
Li 2 [(Lii - X T X )N]; rather, MAE is dominated by single-ion 
anisotropy from impurity T atoms, especially for lower T 
concentration, where T-T atoms become well separated. 
In this work, we investigate the magnetic anisotropy with 
different T elements based on second-order perturbation 
theory by using a Green’s function method. Lorentzians 


1 



are used to represent local impurity densities of states 
and calculate the MAE as a continuous function of band 
filling. First-principles calculations of MAE are also per¬ 
formed to compare with our analytical modeling. 

The present paper is organized in the following way. In 
Sec. |I!| we overview the general formalism of the single¬ 
ion anisotropy[17] 18] with Green’s functions and second- 
order perturbation approach pl9l - f24] . Analytical model¬ 
ing and calculational details are discussed. In Sec. m we 
discuss the scalar-relativistic electronic structure of these 
systems. The band-filling effect on MAE in Li 2 [(Lii_ x T x 
)N], with T=Mn, Fe, Co, and Ni, is examined within 
our analytical model and results are compared with first- 
principles DFT calculations. The results are summarized 
in Sec. IIV1 


II. THEORY AND COMPUTATIONAL DETAILS 

A. Perturbation theory of the magnetocrystalline 
anisotropy and orbital moment 

Perturbation theory allows us to calculate magnetic 
anisotropy directly from the unperturbed band structure. 
Orbital moment, SOC energy, and their anisotropies can 
be written in terms of the susceptibility. [7, 17, 21] [23] The 
relativistic energy E r due to the spin-orbit interaction 
AE so =^L-S can be written as 

1 C Ef rlF 

E r = -~ —Z(Tr[G(E)AV so }) (1) 
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where G(E) is the full Green’s function, which includes 
SOC and can be constructed from the non-perturbed 
Green’s function Go- Using second-order perturbation 
theory (here we consider only systems with a uniaxial 
geometry), the relativistic energy can be written as 


is defined as 

Xmm'( E F) = Xm?m( E F) = [ ~9{3m3m'}i ( 3 ) 

where is the unperturbed on-site Green’s func¬ 
tion. Because we only consider the on-site contribution of 
MAE, only the on-site Green’s function or local suscep¬ 
tibility is needed to investigate MAE. We further assume 
that on-site Green’s functions diagonalize in real har¬ 
monic space. The angular dependence and band struc¬ 
ture dependence of relativistic energy E r are decoupled. 
In the following, we assume that MAE is dominated by 
a particular site i, and consider only its contribution. 

When the spin-quantization axis is along the 001 direc¬ 
tion, the spin-parallel (longitudinal) components of SO 
interaction l z couple orbitals with the same \m\ quantum 
number (m=-m'), while the spin-flip (transverse) ones Z± 
couple orbitals with different \rrt\ numbers (|m| = |m|=bl). 
Hereafter, we refer to those two types of coupling as 
intra-hn| and inter- \m\ types, respectively. According to 
Eq. pi) and absorbing the site index i, the relativistic 
energy can be written as 

e ooi = -j E E (^'C' + ^-o (4) 

CT = ±1 771,771/ 

Positive-definite coefficients A and B are just the spin- 
parallel and spin-flip parts of the \L-S\ 2 matrix elements. 
They can be written as 
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Green’s functions are represented in a basis of or- 
thonormalized atomic functions |i,m,cr), and i labels 
atomic sites, m subbands (in cubic harmonics), and a 
the spin. The local susceptibility Xmm^ characterizing 
the transition between two subbands |ra, a) and |m',cr'), 


A and B correspond to intra-|m| and inter- \m\ transi¬ 
tions, respectively. An interesting property of the coeffi¬ 
cient matrices is 


Brryim/ 


(7) 


For an arbitrary spin orientation other than the 001 
direction, one can either obtain the relativistic energy E r 
by rotating Go [7 or U so ^, 26i in spin subspace. Here we 
use the latter approach and the relativistic energy with 
spin being along the 110 direction can be written as 
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Notice that spin-parallel coefficients in Eq. ® are ex¬ 
actly half of the spin-flip coefficients in Eq. Q. If the 
susceptibility matrix x is relatively homogeneous with re¬ 
spect to spin, then according to Eqs. Q, ([7]), and (|8|, we 
should expect the spin-flip components of the relativistic 
energy E r to be about twice as large as the spin-parallel 
components 27 . This is true for the weakly magnetic 
atoms in different compounds. 


Let us define the orbital moment anisotropy (OMA) 
and MAE, respectively, as Kl = (L z )ooi — (L 2 )no and 
K = E[ 10 -ES OV In this definition, a positive K indicates 
that the system has a uniaxial anisotropy. If Kl is also 
positive, then the system has a larger orbital magnetic 
moment along the easy axis. Using Eq. @ and Eq. (©, 
the MAE K can be written as 


K = i El- 4 " 


n,' B r) 


')(x 


rrim' 


ft 

mm' 


mm 


ft 

mm 1 


/)■ 


(9) 


MAE is resolved into allowed transitions between all 
pairs of orbitals |ra, cr)e^|m / , cr), corresponding to the 
Xmm' f erms - Since A and B are positive definite, the 
coefficient of Xmm' * s positive when ( m=—m ' and g=g') 
or (|m| = |m'| d= 1 and g——g'\ and is negative when 
(■ m=—m' and cr=—cr') or (|m| = |m'| =b 1 and g—g'\ In 
general, the local susceptibility Xmm' * s a ^ so positive def¬ 
inite; hence we have the following simple selection rule 
for MAE: For intra-|ra| orbital pairs, transitions between 
same (different) spin channels promote easy-axis (easy- 
plane) anisotropy; for inter- \m\ pairs, the sign is the other 
way around, i.e., transitions between same (different) 
spin channels promote easy-plane (easy-axis) anisotropy. 
This simple rule is illustrated in Fig. [lj 





FIG. 1: (Color online) Illustration of the dependence of the 
easy-axis direction on the orbital quantum numbers ( m,m') 
and the spin quantum numbers (cr, a') of two subbands. Con¬ 
figurations (a) and (d) favor uniaxial anisotropy, while (b) 
and (c) favor easy-plane anisotropy. The vertical dotted line 
corresponds to the Fermi energy, Ep . The horizontal line 
separates the majority (up) and minority (down) spin chan¬ 
nels. Occupied states with different \m\ numbers are filled 
with different colors. 


Similarly, the OMA Kl can be written as 
Kl = 7^ (Krnrn' ~ ^mm')(Xm m ' — Xmm') 

m,m' 

Hence, OMA originates from the difference between ^ 
and 44- components of each pair susceptibility, while MAE 
originates from the difference between the spin-parallel 
and spin-flip components. If we sum over contributions 
from all the spin components from each pair of orbitals 
( 771 , 777 /) and define 
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then Eqs. and ( [TO] ) can be written as 


K r\ 'y V {Amm' 


(13) 


d^L ~ 2 (Amm' ^mm')Xmm' (14) 

m,m' 

Obviously, the correlation between OMA and MAE [28] 
only happens when the susceptibility is dominated only 
by one of the spin-parallel components. If it is dominated 
by then the system has a smaller orbital moment 
along the easy axis [27]. If it is dominated by y^, then 
the system has a larger orbital moment along the easy 
axis and we have K=^Kp- 

Equation ^ is useful to explain the MAE in two ex¬ 
treme cases, (i) Nonmagnetic limit: Since the orbitals 
are spin independent, we have = Xmm' = Xmm' = 

Xmm’ ■ Xmm’ vanishes for every pair of subbands mm 1 be- 
cause the spin-parallel components cancel out the spin- 
flip ones, (ii) Zero crystal-field limit: Since orbitals are 
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degenerate, S mm /(Anm' — E> mm ')Xra^' in Eq. (N) van- Using the expressions of coefficients in Eqs. © and 
ishes for each of the four spin components era'. Thus the ©, for a d-orbital system, Eq. © can be written as 
total anisotropy vanishes as in a free atom. 


4 3 1 

£ 2 K = 4 X e - 2,2 + X- 1,1 - 2 (x-l ,0 + Xo,l) - 2 (X- 2,-1 + X- 2,1 + X— 1,2 + Xl, 2 ) 


r 


(15) 


where the ordering of the states is \-2)=d xy , |-1 )=d yz , 
\0)=d z 2, |1 )=d xz , and \2)=d x 2_ z 2. Different point-group 
symmetry results in different orbital degeneracy on site 
i. By summing up the coefficients of equivalent orbital 
pairs, Eq. (15) can be simplified. 


For tetragonal, square planar, or square pyramidal ge¬ 
ometries, one pair of orbitals ( d XZl d yz ) is degenerate. 
Equation (15) can be written as 


= 4xl 22 + Xu 


Xl2 3%0 -l X-2,1- 


(16) 


We recover Eq. (13) in Ref. (21 j. 

For linear, trigonal, petagonal bipyramidal, and square 
antiprismatic geometries, besides ( d xz ,d yz ) orbitals, 


(i d x 2_ y 2 ,d xy ) orbitals are also degenerate. Equation (16) 
can be further simplified as 


~^K = 4 X 22 + Xu - 3xoi - 2 X! 2 (17) 


On the other hand, for tetrahedral and octahedral ge¬ 
ometries, five d orbitals split into two groups E g and n g , 
namely, (d z 2 , d x 2 _ y 2 ) and ( d xy , d xz ). One can easily 


show that the right side of Eq. (15) vanishes as expected 
for cubic geometry. 

Similarly, with the coefficient matrices and orbital de¬ 
generacy, one easily recovers the formulas for the orbital 
moment in the tetragonal system as in Ref. m or A1 and 
A2 as in Ref. [7]. 


B. band-filling effect on MAE in a two-level model 

As shown in Eq. (|9|, the MAE and OMA can be re¬ 
solved into contributions from allowed transitions be¬ 
tween all pairs of orbitals. The sign and weight of the 
contribution are determined by coefficients A m?m / and 
Bm,m' , which only depend on the orbital characters of the 
corresponding orbital pairs. On the other hand, Xmm'i 
or its four components Xmm' > are determined by the elec¬ 
tronic structure, namely, the Fermi level (electron occu¬ 
pancy or band filling), band width, crystal-field splitting, 
and spin splitting. Here we investigate the band-filling 
effect on the MAE contribution from a single pair of or¬ 
bitals. For each orbital pair mm', there are four spin 
components: two spin-parallel (tt an d ii) terms and two 
spin-flip terms (ti and . As assumed in the Anderson 





FIG. 2: (Color online) (a) Schematic Lorentzian-shape densi¬ 
ties of states for subbands mn and m !. (b) Xmm' an d its four 
spin components as functions of Fermi energy. The ampli¬ 
tudes of Xmm' with (c) the maximum at e^ ,3 ' > and (d) the 
minimum at as functions of spin splitting As and crystal- 
field splitting Ac. 


model, Lorentzians are used to represent the local densi¬ 
ties of state (LDOS) in our analytical model to illustrate 
the electronic structure dependence of Xmm' an d MAE. 
Similarly, Ebert et al. m used Lorentzians DOS to an¬ 
alytically investigate the orbital magnetic moment and 
relate it to the impurity density of states at the Fermi 
level. For simplicity, we use the same width for every 
Lorentzian orbital, and the on-site Green’s function for 
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subband | m) in one spin channel a is given by 
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E - s<L + iw 


further assume that the two subbands have the same spin 
splitting, 8^ n —8^ n =8^ n , — 8°^, = As, or, equivalently, have 
( 18 ) the same crystal-field splitting, 8 a rn — 8^=8^ —8^, = Ac, 
in the two spin channels. 


where 8 ^ is the band center and w is the half width. The 
corresponding LDOS for subbands | m) and \m') in two 
spin channels are shown in Fig. |2|a). For simplicity, we 


According to Eq. ([3]), the pairwise local susceptibility 
for orbitals |ra,cr) and | m',cr') can be written as 
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X™ m '(Ep) is a positive-definite function for any Ep 
and reaches the maximum at E F =(8^ n + s^,)/ 2. The 
maximum value increases as the two band centers ap¬ 
proach each other until becoming degenerate, because 
the energies required to transfer electrons from occupied 
states to the unoccupied states become smaller. Band 
narrowing increases Xmm' quickly (nearly 1 /w) until it 
reaches the atomic limit. When the bandwidth becomes 
comparable to or smaller than the SOC constant, SOC 
can lift the orbital degeneracy and shift two states, i.e., 
one above and the other below the Fermi level Ep com¬ 
pletely. On the other hand, if the Fermi level sits between 
two well-separated narrow subbands and bandwidth is 
small compared to the distance between the Fermi level 


and the two band centers, w<^Ep — 8 ^ and w<^8 


<j — Ep 


according to Eq. (|19J), then Xmm f= ^/( £ m'~ £ m) does not 
depend on the Fermi energy. 

Using Eqs. © and ( [T9] ) , the dependencies of Xmm' 
and its four spin components on the Fermi energy Ep 
are shown in Figknb). There is one minimum at 8^) and 
two maxima at £p ^, with 


M) 


Si + 82 + As 


+ l —— ^(Ac) 2 + (As) 2 + 4 u> 2 

( 20 ) 

The two maximum peaks originate from the two spin- 
parallel terms an d Xmm' 5 while the minimum orig¬ 

inates from the spin-flip terms — (Xmm' A Xmm')' I n 
Eq. (20), each spin component Xmm' has its maximum 
amplitude when the Fermi level is around the middle 
of the corresponding two band centers. The two spin- 

flip components have their maximum values at the same 
( 2 ) 

Fermi level 8 F because we assume that the two orbitals 
have the same spin splittings. Contributions from the two 
spin-flip components become identical when two states 
| m) and \m') are degenerate. 

As shown in Eqs. § and ( fl3] ) , the MAE coefficients 
for intra-|m| ( A ) and inter- |rapter ms (■ -B ) have different 
signs. To have a large uniaxial anisotropy, the Fermi level 


should be around the 8 F ^ or for intra-|ra| orbital pairs 
and effi for inter- \rrt\ orbital pairs. Two orbitals can ac¬ 
commodate four electrons in two spin channels, and e^) 
roughly corresponds to band filling of one, two, and three 
electrons with i= 1, 2, and 3, respectively. Figures [2|c) 
andj^d) shows the maximum amplitude of Xij{^F=8y) 

as functions of crystal splitting Ac and spin splitting As. 
(1 3) 

For Ep—8 F ’ , it requires Ac=0 to align the two sub¬ 
bands in the same spin channel (two subbands becomes 
degenerate). For Ep=8 F \ it requires As=±Ac to align 
the two subbands in different spin channels. 


C. Crystal structures 

Li 2 (Lii_ ;E T c )N crystallizes in the (a-LUN structure 
type, which is hexagonal and with space group P6/mmm 
(no. 191). The unit cell of a-L^N contains one formula 
unit. There are two crystallographically inequivalent sets 
of Li atoms, Li/ (16) and Li// (2c), with 6/mmm and 
—6m2 point-group symmetries, respectively. The Li/ 
atoms are sandwiched between two N atoms and form 
a linear -Li/-N- chain along the axial direction, while 
Li 11 sites have twofold multiplicities and form coplanar 
hexagons which are centered at -Li/-N- chains and paral¬ 
lel to the basal plane. Li// is more close packed in lateral 
directions and 3d atoms randomly occupy Li/ sites. We 
carried out DFT calculations for small doping concen¬ 
tration with x=0.166 and found that all T elements with 
T=Mn, Fe, Co, and Ni indeed prefer to occupy Li/ sites. 
To calculate the electronic structure and MAE, we use a 
supercell which corresponds to a '/3x\/3x2 superstruc¬ 
ture of the original a-L^N unit cell. Details of the super¬ 
cell construction can be found in Ref. [12 . For x=0.5, as 
shown in Fig. [3j there are three T atoms in the 24-atom 
supercell with one on the la site and the other two on 
the 2d sites. Both T\ a and T^d sites are derived from the 
16 site in the original a-L^N. They have a linear geome¬ 
try and a strong hybridization with neighboring N atoms 
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LMTO (FP-LMTO) method to check our calculational 
results. For the MAE calculation, the /c-point integration 
was performed using a modified tetrahedron method with 
Blochl corrections, with 16 3 ^-points in the first Brillouin 
zone of the 24-atom unit cell. By evaluating the SOC ma¬ 
trix elements (Vso) and its anisotropy[27, we resolve the 
anisotropy of orbital moment and MAE into sites, spins, 
and orbital pairs. The correlation effects are also consid¬ 
ered by using the local-density approximation (LDA)+U 
method. Here we choose the fully localized limit imple¬ 
mentations of the double counting introduced by Liecht¬ 
enstein et al.[ 35. considering it is more appropriate for 
materials with electrons localized on specific orbitals. 


III. RESULTS AND DISCUSSIONS 
A. Electronic structures 


FIG. 3: (Color online) Schematic representation of the su¬ 
percell used in the DFT calculation for Li 2 [(Lii_ £C T x )N] with 
£=0.5. Both Tia and T^d sites are derived from the Li/ (15) 
site in the original o-LLN structure, while other Li atoms, 
which form coplanar hexagons, correspond to Li// (2c) sites 
in the original o-LLN structure. 


along the axial direction. T\ a have six Li neighbors, while 
T^d have three T^d and three Li neighbors in the T-Li 
plane. This structure (denoted as hex 2 in Ref. [12]) is 
of particular interest because two types of T sites, T\ a 
and possess very different local surroundings and 
represent different local impurity concentrations. Along 
the in-plane direction, T-T distances are rather large, es¬ 
pecially for the la site. Since the T\ a site represents a 
relatively low impurity concentration and dominates the 
uniaxial MAE for T=Fe, most of the results in this work 
are focused on the T\ a site in the hex2 supercell. We 
also consider other concentrations such as x=0.16 and 
x=0.33. 


D. DFT calculational details 

We carried out first principles DFT calculations using 
the Vienna ab initio simulation package (VASP) [29l [30] 
and a variant of the full-potential linear muffin-tin orbital 
(LMTO) method[3T]. We fully relaxed the atomic posi¬ 
tions and lattice parameters, while preserving the sym¬ 
metry using VASP. The nuclei and core electrons were 
described by the projector augmented-wave potential(32] 
and the wave functions of valence electrons were ex¬ 
panded in a plane-wave basis set with a cutoff energy 
of 520 eV. For relaxation, the generalized gradient ap¬ 
proximation of Perdew, Burke, and Ernzerhof was used 
for the correlation and exchange potentials. The spin- 
orbit coupling is included using the second-variation 
procedure [331 IS]. We also calculated the MAE by car¬ 
rying out all-electron calculations using the full-potential 


Without considering SOC, the axial crystal field on 
both T\ a and T^d sites splits five 3d orbitals into three 
groups: degenerate ( d xy , d x 2 _ y 2 ) states, degenerate ( d yz , 
d xz ) states, and d z 2 state. Equivalently, they can be la¬ 
beled as m=z L2, ra=±l, and m=0 using cubic harmonics. 

The scalar-relativistic partial densities of states 
(PDOS) projected on the site are shown in Fig. |4| 
For T=Fe, the PDOS obtained is very similar to what 
was previously reported [12. The Fe 3d shell has seven 
electrons and the majority spin channels of d orbitals are 
fully occupied with five electrons. 

The Fe d z 2 states hybridize with p z states of N atoms 
along the axial direction and mix with on-site 4s states, 
which causes the d z 2 orbital to be lower in energy than 
the other d orbitals. [12] The d z 2 states spread out and 
lie below the Fermi level and accommodates one electron 
in the minority spin channel. The last electron occupies 
half of the degenerate (d xy , d x 2 _ y 2 ) states in the minority 
spin channel. These states have a very narrow bandwidth 
and cross the Fermi level. 

The linear geometry minimizes the in-plane hybridiza¬ 
tion between the T 3d orbitals and the neighboring 
atoms, making them atomic like and resulting in nar¬ 
rower bands. The T^d site shows a similar PDOS as the 
T la site; however, the in-plane hybridization with other 
T 2 d sites results in a much broader bandwidth than the 
la sites. 

For other T elements, the DOS peaks are well sep¬ 
arated as in T=Fe. The minority spin channel clearly 
shows a different band-filling pattern with different T el¬ 
ements. The deviation from the rigid-band model is also 
obvious. Spin splitting decreases from Mn to Ni, while 
the crystal-field splitting values (the energy difference be¬ 
tween ra==bl and ra==b2 states) are larger for T=Mn and 
Fe than for T=Co and Ni. 

Figure [5][a) shows the schematic Fe PDOS, and how 
the Fermi level changes with different T in a rigid-band 
approximation (RBA). Different T elements correspond 
to different integer number of 3d electrons. Since each 
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FIG. 4: (Color online) Partial densities of states projected 
on the 3 d states of the T\ a site in the hex 2 structure in 
Li 2 [(Lii_ x T x )N], where x=0.5 and T is (a) Mn, (b) Fe, (c) Co, 
and (d) Ni. The vertical dotted line corresponds to the Fermi 
energy, Ef . The horizontal dotted line separates the major¬ 
ity (up) and minority (down) spin channels. Calculation is 
within LDA, without spin-orbit coupling included. 


degenerate state pair can accommodate two electrons in 
one spin channel, the Fermi level either intersects the 
degenerate peaks or sits in the middle of two peaks. 


B. MAE in Li2[(Lii_ a; T £C )N] with T= Fe 

MAE in Li 2 [(Lii_ ;E T E )N] with T=Mn, Fe, Co, and Ni 
and x = 0.5 are calculated in DFT and summarized in 


TABLE I: Lattice constants, total and site-resolved MAE 
Li 2 [(Lio. 5 To. 5 )N] with T=Mn, Fe, Co and Ni. The MAE val¬ 
ues for T 2 d site are in unit of meV/ atom, and there are two 
T^d atoms in the supercell. 


T 

Lattice parameters 

K{meV) 


a(a.u.) 

c/a 

cell Tia 

T2d 

others 

Mn 

12.143 

1.202 

-1.14 -0.35 

-0.38 

-0.03 

Fe 

12.091 

1.183 

20.83 14.77 

3.09 

-0.12 

Co 

12.144 

1.154 

-3.69 -0.89 

-1.32 

-0.15 

Ni 

12.113 

1.156 

2.52 1.71 

0.37 

0.06 


Table|IJ The system has uniaxial anisotropy with T=Fe or 
Ni and easy-plane anisotropy with T=Mn or Co. MAE is 
dominated by the contributions from the la site for T=Fe 
or Ni. Results are in qualitative agreement with previous 
calculations. pTHl3j The extraordinary MAE for T =Fe 
originates from the unique band structure in this system. 
Because the well-isolated Fe atoms, such as the Fei a site 
in the hex2 supercell, provide the major contribution to 
the uniaxial anisotropy, we focus on the Fei a site. 

As shown in Fig. [l] the sign of the MAE contribu¬ 
tion from transitions between a pair of subbands |ra, a) 
and | m',a') is determined by the spin and orbital char¬ 
acter of the involved orbitals. Because the d z 2 orbital is 
spread out relatively further below the Fermi level and 
contributes negligibly to the MAE, we only consider the 
transitions between subbands with m= —2, —1, 1, and 2. 
Intra-|ra| transitions |1)<R>| — 1) and |2)«a| — 2) promote 
easy-axis anisotropy when they are within the same spin 
channel, and easy-plane anisotropy when between dif¬ 
ferent spin channels. For inter- \m\ transitions, it is the 
other way around. Transition | db 1)<B>| =b 2) promotes 
easy-plane anisotropy when it is within the same spin 
channel and easy-axis anisotropy when between different 
spin channels. The signs and coefficients of the MAE 
contributions from different orbital pair transitions are 
indicated in Fig. |5|a). Transitions contribute to MAE 
only when they cross the Fermi level. The amplitude 
of MAE depends on the orbital characters and also the 
energy difference between the two band centers. When 
the Fermi level intersects the narrow degenerate states, 
the transition energy required to excite an electron across 
the Fermi level is very small (between 0 and bandwidth), 
making the MAE contribution from this pair of orbitals 
very large. On the other hand, when the Fermi level 
is between two well-separated DOS peaks, the required 
transition energy is much larger so the amplitude is much 
smaller. 

To elucidate the orbital contributions from the Fei a 
site to the MAE in Li 2 [(Li 0 . 5 Feo. 5 )N], we approximate 
the densities of states of | d= 1) (d xz ,d yz ) and | d= 2) 
( d xy ,d x 2_ y 2 ) subbands with two Lorentzian functions. 
Crystal-field splitting Ac=e| m | = i — e| m |= 2 =1.8eV, spin 
splitting As=2.4 eV, and half width w = 0.06 eV are 
used to represent the DFT-calculated PDOS, as shown in 
Fig. g The PDOS used in our model is shown in Fig. Ha) 
and the MAE contribution from the la site and its de- 


7 

















































E(eV) 



^F(eV) 



FIG. 5: (Color online) (a) Schematic partial densities of states 
projected on the 3 d states of Fei a sites. Orbital transitions 
and the sign of their contributions to the MAE are also shown. 
Solid line indicates positive contribution (easy axis) and the 
dashed line indicates negative contribution (easy plane) to 
the easy-axis anisotropy, (b) Scaled MAE 4 K/Z 2 from T\ a 
site and its decomposition into orbital susceptibilities as func¬ 
tions of band Filing, (c) Magnetic anisotropy energy K from 
T\ a site as a function of T. Different sets of electronic struc¬ 
ture parameters As, Ac, and w are used to represent the 
DFT PDOS on T\ a sites in Li 2 [(Lio. 5 To.s)N] for different T 
elements. 


composition into orbital pair transitions as functions of 
the Fermi energy are shown in Fig.[5|b). With T=Fe, the 
Fermi level intersects the | ±2,^) states, which results in 
a large uniaxial anisotropy. Using Eq. ( p~7| ), Fei a has a 
MAE contribution which is of the order of 15 meV/Fe. 
As shown in Fig.[5|b), for T— Fe, nearly all MAE contri¬ 
butions are from the transitions |2,^)<e>| — 2,4,), in other 
words, between d x 2_ y 2 and d xy orbitals in the minority 
spin channel. 

To compare with the above analytical modeling, MAE 
calculations were carried out in both VASP and all¬ 
electron FP-LMTO. The difference of MAE values using 
two methods is less than 5% for T=Fe. To decompose 
the MAE, we evaluate the SOC matrix element ( V so ) and 
its anisotropy iF((U so )), which can be easily decomposed 
into sites, spins, and orbital pairs[27:. We found that 
K « K((V so ))/2 for all T compounds, which suggests 
that second-order perturbation theory is a good approx¬ 
imation. As shown in Table [TJ for T=Fe, the total MAE 
is 20.8 raeV (per 24-atom cell) and MAE contributions 
from la and 2 d sites are 14.77 and 3.09raeV/Fe, respec¬ 
tively. The contributions from Li and N atoms are nearly 
zero as expected. Thus, the impurity Fe (especially Fei a ) 
atoms are essentially the only MAE providers. By fur¬ 
ther investigating the matrix element of SOC on the la 
site, we found that nearly all the MAE contributions 
came from intra-|m| transitions of |2,|)o| — 2,|). As 
shown in Table |ll[ the 4^22 term (dominated by X 22 f° r 
T=Fe ) contributes 15.1 meV /Fe and the xh term has a 
much smaller negative value of -0.42 meV /Fe, while other 
terms are negligible. Hence, DFT results agree with our 
model very well. 

With magnetization along the c direction, the SOC can 
lift the orbital degeneracy and shift two narrow bands 
m = ±2, one below and the other above the Fermi level 
completely, with orbital quantum number m c = ±2, re¬ 
spectively, where m c is the orbital quantum number in 
the complex spherical harmonics. As a result, the density 
of states at the Fermi level becomes very small. Indeed, 
experiments [15] found this system to be an insulator for 
T=Fe. It had been shown that [HHI3ES1 the correlation 
effect further enhances the separation between occupied 
and unoccupied states. Using the LDA+U method, we 
also found that correlation can enhance the orbital mo¬ 
ment when the spin is along the axial direction. 

Fe concentration and site disordering can significantly 
affect the MAE. As we have shown, the Fe 2 d sites, which 
represent a high-doping concentration, have much lower 
anisotropy than the Fei a sites, which represent a lower- 
doping concentration. By replacing the Fe 2 d sites back 
with Li atoms in the hex 2 supercell, we calculated the 
MAE with a smaller concentration x=0.166 and found 
that MAE increase to 22 meV/Fe, which is in very good 
agreement with previous calculations. m- An interest¬ 
ing concentration is x=0.33. If only one of two 2d sites is 
occupied by Fe in the hex2 supercell, as shown in Fig. [3j 
then this configuration would correspond to x=0.33 and 
the supercell has two well-isolated Fe atoms. The DFT 



























calculation shows high MAE with a value of 20 meV/Fe. 
On the other hand, if the two Fe atoms occupy the 2d 
sites and then are not well separated, the resulting MAE 
is much smaller (2.8 meV/Fe). Even if we assume that 
Fe atoms tend to separate, with a concentration beyond 
x=0.33, it is unavoidable to have Fe atoms neighboring 
each other and the hybridization between them causes 
the MAE (per Fe) to decrease. Furthermore, impurity 
sites are disordered, as found in experiments. At least at 
a higher concentration, many Fe atoms would not have 
the symmetric lateral surroundings as the two Fe sites do 
in the hex 2 supercell we used in the calculations. This 
site disordering may also have an effect on MAE by lower¬ 
ing the point-group symmetry of Fe impurity sites. And 
the m = ±2 states on Fe sites are no longer degenerate, 
which may decrease MAE per Fe. 


C. MAE in Li2[(Lii_ x T x )N] with T=Mn, Co and 
Ni: The band-filling effect 

Figure [5|a) shows how the Fermi level changes with 
different T elements in a simple rigid-band picture. Only 
those transitions across the Fermi level contribute to 
MAE. With T elements other than Fe, the | ±2,4,) 
states become either fully occupied or unoccupied. The 
large uniaxial anisotropy that originated from transition 
12, i)^\ — 2, 1) (term 4x^2) vanishes and other transitions 
becomes important, depending on the position of the 
Fermi level. For T=Ni, the Fermi level intersects the de¬ 
generate | ± 1,|) states. Hence anisotropy contributions 
are dominated by the transitions |1,|)<R>| — 1,|) (term 
xfi)- This transition promotes the uniaxial anisotropy, 
as 4%22 does f° r T=Fe. For T = Co, the Fermi level is 
between | ±2,^) and | ± 1,^) peaks. The transitions of 
| ±2,|) | ±1,|) (term -3%j^) and | ± 1, t) | T 1,1) 

[term -(Xn+xfl)] support easy-plane anisotropy, while 
the transition | ± 2,t)^| ± 1,|) (term 3 X 12 ) promotes 
easy-axis anisotropy. However the two bands involved in 
the last transition are far away from each other and this 
contribution is relatively small. Hence, for T=Co, one 
should expect the system to have easy-plane anisotropy. 
For T=Mn, there are four transitions that contribute to 
the MAE; all of them are between the two spin channels, 
in which two inter- \m\ transitions | ± 1, t)^| ±2, \) (term 
3 X 12 ) and | ±2, t)^H ± 1? i) support easy-axis anisotropy, 
while two other intra-|m| transitions |±1, t)«-H=Fl, -1} and 

I ± 2 , t)«-H T 2 , 4 .) [term - 4 (x 22 +X 22 )] support easy-plane 
anisotropy. The four transitions compete and the sign of 
the total MAE is not obvious and requires a more quan¬ 
titative description. 

The SOC constant £ changes with element. In 
Fig. [^b), we plot the scaled MAE K=K/ 4£ 2 and its 
orbital-resolved components as functions of the Fermi 
level by using parameters of As, Ac, and w for T=Fe. In 
a rigid-band picture, it clearly shows that Ni also has 
a uniaxial anisotropy with contributions coming from 


the xfi term. Since we are using the same half width 
w of LDOS for ra=±l and m=±2 subbands, we have 
Kjsfi ~ \Kfc because of the intra-|m| transitions coef¬ 
ficients m 2 , as shown in Eqs. ^ and Figure [HJc) 
shows the MAE IT as a function of the number of oc¬ 
cupied electrons by using different sets of As, Ac, and 
w parameters to better present DFT-calculated PDOS 
for different T elements, as shown in Fig. [4j The SOC 
constant £ is interpolated by using DFT-calculated £ val¬ 
ues for 3d elements. Since £ decreases with the atomic 
number within a given nl shell, K quickly decreases with 
smaller atomic numbers due to the factor £ 2 . The DFT 
MAE values are also plotted to compare with the mod¬ 
eling MAE function. As shown in Fig. [SJc), with T=Fe 
parameters, the modeling MAE (Fe rigid-band approx¬ 
imation) can already correctly describe the MAE trend 
with different T elements. 

Although the RBA predicts the correct easy-axis di¬ 
rection for T=Ni, the difference between RBA mod¬ 
eling and DFT is rather large. In RBA modeling, 
Kni/K pe^i^Ni/^Fe ) 2 /4 ~0.6, while the DFT value 
(1.71raeF/atom) for T=Ni is about one order of mag¬ 
nitude smaller than for T=Fe. This can be explained 
as follows. First, we use the same band width for all 
DOS peaks in our modeling. In fact, the | ± 1,4-) bands 
are much broader than the | ± 2,|) bands. The easy- 
axis anisotropy contribution from the transition between 
| ± 1, l) states decreases with increasing band width. Sec¬ 
ond, the Ni PDOS deviates from the Fe PDOS more than 
Mn or Co, so RBA is less appropriate for T=Ni. The 
spin splitting As and crystal-field splitting Ac are much 
smaller in Ni than in Fe. This causes the amplitudes of 
the negative contributions from | ± 2,^)e^| ± 1,^) and 
| ± 1, t)^| T1,1) to become larger and decrease the total 
uniaxial anisotropy. As shown in Fig. (5|c), if we use a 
smaller As, smaller Ac, and larger w to better represent 
the Ni PDOS calculated from DFT calculations, then 
much better agreement between model and DFT values 
can be reached. 

For T=Co, the model MAE is about twice the DFT 
value, probably because of the simplified model DOS. 
The orbital-resolved T\ a MAE calculated in DFT are 
summarized in Table [TTJ Overall, there is a qualitative 
agreement between DFT and the analytical model for 
the orbital-resolved MAE values for all T elements. It 
is interesting that with T=Co, the contribution of the 
4X22 term is comparable to that of — 2 %i 2 and Xii in 
DFT, which is not expected in the model. As shown in 
Fig. [4](c) , there is a small portion of unoccupied | ± 2,|) 
states right above the Fermi level in the minority spin 
channel, which makes the 4%22 terms comparable to oth¬ 
ers. However, this electronic structure detail is not con¬ 
sidered in the simplified DOS we use in modeling. If we 
neglect the 4%22 terms in DFT, then a better agreement 
between modeling and DFT can be achieved for T=Co. 

Thus, the contributions from well-separated impurity 
sites with T can be well understood. For T=Mn and 
Co, the easy-plane anisotropy is a result of competition 
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TABLE II: Orbital-resolved MAE from the T\ a site in 
Li 2 [(Lio. 5 To. 5 )N] with T— Mn, Fe, Co, and Ni. 







K (meV) 


Term 

Orbital Transition 

Mn 

Fe 

Co 

Ni 

4X22 

dxy 


d x 2 — 

-0.86 

15.10 

0.71 

-0.03 

Xii 

dyz 


d xz 

-0.22 

-0.42 

-0.78 

3.68 

-2xi2 

dyz )d X Z 


dxy id x 2_ y 2 

0.73 

-0.18 

-0.81 

0.09 

-3xoi 

d z 2 


dyz ,dxz 

0.03 

0.08 

-0.01 

-0.25 


between different transitions, instead of being dominated 
by the intra-|m| transition, which strongly depends on 
the bandwidth of the degenerate | d= m) states that are 
intersected by the Fermi level. As a result, the band¬ 
narrowing effect on MAE is not as strong as for T =Fe or 
Ni. As shown in Table [I] the contributions from 2d sites 
are comparable or even larger than la sites for T=Mn 
and Co. 


IV. SUMMARY AND CONCLUSION 

Based on second-order perturbation theory, MAE is 
resolved into contributions from different pairs of or¬ 
bital transitions, more precisely, the difference between 
spin-parallel and spin-flip components of the orbital sus¬ 
ceptibilities of the corresponding orbital pair. In the 
Li 2 [(Lii_ a ,T a ,)N] systems, with T=Mn, Fe, Co, and Ni, 
the linear geometry of the T sites minimizes the in-plane 
hybridization and results in atomic like orbitals around 
the Fermi level for all T elements. The MAE oscillates 
with the atomic number from T=Mn to T=Ni, which is 
a result of the competition between contributions from 


all allowed orbital transitions. As the Fermi level evolves 
with T, different orbital pair transitions dominate the 
contribution to MAE. For T=Fe and T=Ni, the intra-|ra| 
transitions within the minority spin channel dominate 
the MAE contribution and result in a uniaxial anisotropy. 
For T=Mn and Co, the easy-plane anisotropy is a result 
of the competition between contributions from several 
transitions with different signs. Using Lorentzian den¬ 
sity of states, we investigate the band-filling effect on 
MAE in an analytical model based on a Green’s function 
technique. We show the MAE as a continuous function 
of atomic number. This analytical model can already 
describe the correct trend of the MAE obtained using 
DFT, by just using a simple rigid Fe band picture. If we 
take into account the deviation from the rigid Fe band 
model and some details of DFT electronic structure, an 
even better agreement between the model and DFT can 
be found. To further validate our modeling analysis, we 
also calculate the orbital-resolved MAE by evaluating the 
SOC matrix element in DFT. Overall, Li 2 [(Lii_ iC T E )N], 
with T=Mn, Fe, Co, and Ni, is a unique system which 
clearly shows the band-filling effect on MAE and the na¬ 
ture of this effect can be understood in a very simple 
model. 
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